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A  COMPUTATIONAL  STUDY  OF  THE  CHEMICAL  KINETICS 
OF  HYDROGEN  COMBUSTION 


I.  INTRODUCTION 

A  major  difficulty  in  the  study  of  reactive  flows  is  the  investigation  of  the  influence  of  various 
contributing  processes  in  a  system's  behavior.  Because  of  this,  it  is  often  desirable  to  examine  the 
component  processes  separately  before  attempting  to  understand  their  interaction.  In  combustion  sys¬ 
tems,  for  example,  it  is  important  to  investigate  the  chemical  and  gasdynamic  processes  individually 
before  considering  the  fully-coupled,  nonlinear  system.  However,  the  strong  coupling  between  these 
processes  sometimes  poses  an  added  problem  to  their  separate  and  independent  study. 

fn  addition  to  these  problems,  there  are  other  difficulties  that  are  associated  more  explicitly  with 
the  study  of  the  individual  processes.  One  important  problem  that  arises  in  the  study  of  the  chemistry 
of  combustion  systems  is  the  problem  of  accurately  determining  the  reaction  mechanism  that  describes 
the  system  over  the  range  of  temperatures,  pressures  and  stoichiometries  encountered  during  combus¬ 
tion.  Such  a  reaction  mechanism  is  usually  composed  of  a  set  of  elementary  or  primary  reactions  that 
involve  a  number  of  intermediate  species  that  are  produced  by  initial  reactants  in  their  conversion  to 
products.  However,  even  in  the  simplest  chemical  kinetic  schemes,  it  is  often  difficult  or  nearly  impos¬ 
sible  to  determine  all  of  the  elementary  reactions,  their  rates  and  intermediate  species.  Experiments 
can  determine  the  products  of  combustion  and  can  sometimes  detect  significant  metastable  intermedi¬ 
ates  and  define  their  related  reaction  cross  sections,  but  inference  and  supposition  are  always  involved 
in  the  assignment  of  a  reaction  mechanism. 
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Our  efforts  in  modeling  combustion  in  hydrogen-containing  gas  mixtures  have  led  us  to  consider 
the  accuracy  with  which  a  proposed  H2  -  0:  elementary  reaction  set  describes  the  chemistry  of  this  sys¬ 
tem.  These  considerations  are  important  because  hydrogen  itself  is  an  important  fuel  and  because  its 
oxidation  mechanism  plays  a  fundamental  role  in  the  combustion  of  hydrocarbons.  An  accurate 
description  of  hydrogen  combustion  is  also  useful  because  it  will  help  in  the  evaluation  of  the  likelihood 
of  deflagration  and  explosion  hazards  in  industrial  facilities  such  as  nuclear  reactors.  Therefore,  the 
purpose  of  this  paper  is  to:  1)  present  a  homogeneous,  gas-phase  reaction  mechanism  descriptive  of  the 
Hj  — 02  combustion  system,  2)  describe  the  computational  kinetic  work  that  we  have  done  to  test  the 
accuracy  of  the  mechanism,  and  3)  summarize  our  conclusions  regarding  the  general  applicability  of 
this  mechanism. 

In  the  work  presented  here,  we  have  assumed  that  the  ,  ✓iduai  states  of  the  reacting  molecules 
are  equilibrated  and  that  local  thermodynamic  equilibrium  exists.  This  simplifies  the  problem  and 
allows  us  to  use  the  principle  of  thermodynamic  consistency  between  forward  and  reverse  reactions  to 
compute  poorly  known  reaction  rates  via  the  equation: 

f^eq  " "  k  forward!  ^ reverse  •  Ilf 

Here  Kcq  is  the  thermodynamic  equilibrium  constant  and  k,nrnard  and  k,eimr  are  the  reaction  rate 
coefficients  for  the  associated  forward  and  reverse  reactions. 

In  order  to  test  the  proposed  reaction  mechanism,  a  series  of  detailed  chemical  kinetic  calculations 
have  been  performed.  The  computational  code,  referred  to  as  CHEMOD,  has  been  described  in  detail 
by  Flanigan  and  Young  [1].  Briefly,  the  code  consists  of  three  principal  parts:  1)  an  I/O  and  logical 
control  program,  2)  a  syntactical  reaction  transcriptor  (SRT)  and  differential  function  generator,  and  3) 
a  vectorized  version  of  the  asymptotic  integration  method  developed  by  Young  [21.  The  SRT  was 
developed  by  Young  explicitly  for  the  CHEMOD  code.  The  code  computes  total  molar  concentration, 
total  mass  density,  total  pressure,  internal  energy,  and  individual  species  concentrations  as  a  function  of 
time  under  isothermal  or  adiabatic  constant-volume,  ideal-gas  reaction  conditions. 
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The  results  obtained  from  these  calculations  have  been  compared  with  available  measurements  of 
reaction  induction  times,  second  explosion  limits,  the  temporal  behavior  of  reaction  species  and  reac¬ 
tion  equilibrium  parameters.  These  comparisons  have  been  critically  evaluated  with  regard  to  the 
experimental  conditions  of  the  measurements.  Specifically,  contributions  by  gasdynamic  effects  have 
been  noted. 

II.  THE  H:  -  O,  REACTION  MECHANISM 

Ever  since  the  early  work  of  Hinshelwood,  Sagulin  and  Semenov  [3-7),  the  H2-O2  reaction  sys¬ 
tem  has  been  the  subject  of  continuous  and  intensive  research  [8-13],  In  this  section,  we  summarize 
some  of  the  important  features  of  this  system  in  order  to  provide  a  background  for  the  material  that 
follows.  In  addition,  we  discuss  some  of  the  details  pertinent  to  the  mechanism  that  we  have  used  to 
describe  the  H:-0:  reaction  system. 

The  gas-phase  reaction  between  H2  and  O;  is  a  complex  reaction  system  composed  of  a  relatively 
large  number  of  elementary  reactions.  The  system  consists  of  a  chain  reaction  sequence  that  branches 
between  the  atomic  and  free  radical  chain  carriers  H.  O.  HO  and  HO:.  At  low  temperatures.  HO:  is 
only  moderately  reactive  and  is  important  in  chain  termination.  At  temperatures  above  —  1500°K.  HO: 
is  more  reactive  and  becomes  an  effective  chain  carrier. 

A  principal  feature  of  the  H2-0:  reaction  system  is  its  demonstration  of  explosion  limits  that 
vary  as  a  function  of  pressure  and  temperature.  Information  about  these  limits  was  initially  deduced 
from  observations  made  in  static  cell  experiments  [14-22].  These  studies  showed  that  cell  dimensions 
and  wall  surface  materials  are  important  factors  and  led  to  the  conclusion  that  heterogeneous  reactions 
make  important  contributions  in  the  H2  -  0-  system.  This  is  especially  true  at  the  first  explosion  limit. 
The  second  explosion  limit  is  also  affected,  but  to  a  much  lesser  extent. 


3 


BURKS  AND  ORAN 


Figure  1  shows  the  explosion  limit  diagram  as  a  function  of  pressure  and  temperature  for  a 
stoichiometric  mixture  of  Hi  and  02.  This  figure  can  be  qualitatively  understood  in  terms  of  the  rela¬ 
tionship  between  chain  branching  and  chain  termination  first  pointed  out  by  Semenov  [9]  and  Hinshel- 
wood  [8].  They  have  shown  that  whenever  the  rate  of  chain  branching  exceeds  the  rate  of  chain  termi¬ 
nation,  fast  combustion  occurs.  Inversely,  whenever  the  rate  of  chain  termination  exceeds  the  rate  of 
chain  branching,  H2  and  02  react  slowly. 

Thus,  at  initial  pressures  below  the  first  explosion  limit  (where  the  mean  free  path  of  the  chain 
carriers  is  comparable  to  reaction  vessel  dimensions),  chain  termination  by  heterogeneous  wall  reactions 
predominates  and  a  slow  reaction  is  observed.  As  the  pressure  of  the  reaction  mixture  is  increased 
(and  the  mean  free  path  of  the  chain  carriers  is  decreased),  the  rates  of  the  bimolecular  chain  branch¬ 
ing  and  propagation  reactions  are  increased  relative  to  the  heterogeneous  chain  termination  reactions. 
Finally,  chain  branching  exceeds  chain  termination  at  the  first  explosion  limit  and  the  overall  rate  of  the 
reaction  increases  exponentially. 

Between  the  first  and  second  explosion  limits,  in  the  region  of  the  explosion  peninsula,  the  gas 
mixture  explosively  ignites  at  all  pressures.  At  the  second  explosion  limit,  chain  termination  by  ter- 
molecular  reactions  becomes  important.  As  the  pressure  of  the  reaction  mixture  is  increased  to  just 
above  the  second  explosion  limit,  the  rate  of  termolecular  chain  termination  overbalances  the  rate  of 
chain  branching  and  a  slow  reaction  between  H2  and  02  again  predominates.  The  overall  rate  of  reac¬ 
tion  between  H2  and  02  then  decreases  to  a  minimum  and  then  begins  to  slowly  increase  as  the  pres¬ 
sure  is  raised  above  the  second  explosion  limit.  The  overall  reaction  rate  continues  to  increase  with 
increasing  pressure  until  the  rate  of  thermochemical  heat  release  exceeds  the  rate  of  heat  removed.  In 
this  region,  the  reactions  become  self-accelerating,  or  autocatalytic,  and  the  third  explosion  limit 
becomes  defined.  At  pressures  above  the  third  explosion  limit,  it  is  thought  that  all  pressures  ignite  to 
explosion. 
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The  Hi  — Oi  reaction  mechanism  tested  in  this  work  includes  all  likely  homogeneous,  gas-phase 
reactions  among  all  feasible  neutral,  electronic  ground  state  species.  Reactions  were  not  included  that 
involved  mechanistically  difficult  atomic  rearrangements  equivalent  to  the  operation  of  two  or  more 
other  elementary  reactions.  Where  available  data  showed  that  a  reaction  was  not  well-known,  this  reac¬ 
tion  was  omitted.  The  reaction  mechanism  incorporates  termolecular  reactions  presuming  that  the 
major  third  body  constituent  of  these  reactions  is  a  heavy  inert  species  such  as  argon.  In  Table  1,  we 
present  the  mechanism  as  a  series  of  forward  and  reverse  elementary  reactions  and  rate  constants. 

Reaction  rate  constants  were  chosen  from  a  survey  of  available  data  and  are  expressed  in  the 
modified-Arrhenius  form: 

k(T)  -  ATBexp(C/T)  (2) 

where  k  ( T)  is  the  temperature-dependent  reaction  rate  constant,  T  is  temperature  (°K)  and  A,  B  and  C 
are  constants.  For  those  reactions  for  which  rate  constant  information  was  either  sparse  or  not  known 
at  all,  an  individual  rate  constant  was  computed  using  Equation  (1). 

Concentration-based  equilibrium  constants,  AT/,,  were  cast  into  the  same  modified-Arrhenius  form 
as  the  reaction  rate  constants.  A  third-order  matrix  expression  was  constructed  by  selecting  thermo¬ 
dynamic  equilibrium  constant  data  at  three  points  that  spanned  the  300-3000°K  temperature  range  for 
the  following  equation: 

In  K^qJ  -  In  Aeq  +  Beq  In  T{  -  CJT,  (3) 

where  K$q,  and  T,  are  the  equilibrium  constant  and  temperature  of  the  /  th  point  and  Ae<r  Beq  and  Ceq 
are  constants.  The  thermodynamic  data  was  taken  from  the  JANAF  Tables  [231  and  from  Engleman’s 
compilation  of  kinetic  data  [24],  Solution  of  the  matrix  relation  for  the  three  constants— Aeq,  Bcq  and 
Ceq— for  each  forward  and  reverse  reaction  pair  then  provided  for  the  temperature-dependent  descrip¬ 
tion  of  the  equilibrium  constants  as  shown  in  Table  2. 
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Forward  and  reverse  reaction  rate  constants  were  checked  for  thermodynamic  consistency  over  the 
300-3000°K  temperature  range  by  comparison  of  temperature  plots  of  reaction  rate  constants  obtained 
by  experiment  and  by  use  of  Equation  (1).  When  a  rate  constant  was  either  poorly  known  experimen¬ 
tally  or  was  not  consistent  with  its  reverse  rate  and  equilibrium  constant,  the  computed  rate  constant 
was  used  in  lieu  of  the  experimental  rate  constant.  Since  there  are  frequently  more  problems  and  larger 
errors  associated  with  the  measurement  of  rate  constants  at  elevated  temperatures,  most  of  the  com¬ 
puted  rate  constants  are  among  the  endoergic  reactions  listed  as  reverse  reactions  in  Table  1.  Three 
examples  of  reactions  that  were  found  to  be  thermodynamically  inconsistent  are  shown  in  Figures  2-4. 

The  use  of  the  thermodynamic  equilibrium  relationship  to  compute  unknown  rate  constants  is 
based  on  the  assumption  that  the  elementary  reactions  of  the  H2-02  combustion  system  involve 
reactants,  intermediates  and  products  that  can  be  described  by  Boltzmann  thermal  energy  distributions 
This  condition  is  generally  correct  in  the  slow  reaction  zones  of  concentrated  gas  mixtures  and  even  in 
the  explosively  reactive  zones  of  highly  diluted  gas  mixtures.  The  assumption  is  supported  by  the  work 
of  Asaba  et  al.  [25]  and  Belles  and  Lauver  [26]  who  have  shown  that  vibrational  equilibration  of  02  in 
shock-initiated  combustions  only  affects  the  H2-02  reaction  kinetics  by  a  very  small  amount.  Thus, 
since  it  appears  that  large  shock-induced  disturbances  are  satisfactorily  equilibrated  within  the  induction 
time  of  the  reactive  gas  mixture,  we  presume  that  more  uniform  energy  perturbations  due  to  combus¬ 
tion  (exclusive  of  explosion)  are  efficiently  relaxed  so  that  equilibrium  kinetics  can  be  assumed. 

The  modified-Arrhenius  form  of  reaction  rate  constants  that  is  commonly  used  has  been  obtained 
by  empirical  data  fitting.  Although  transition  state  (or  activated  complex)  theory  [27-29]  and  collision 
theory  [30-32]  predict  non-Arrhenius  temperature  behavior,  neither  theory  provides  an  explicit,  analyti¬ 
cal  temperature  dependence  for  the  rate  constant,  k{T). 

Justification  for  the  use  of  the  modified-Arrhenius  rate  constants  has  been  provided  by  a  number 
of  workers  [27-36].  Zellner  [28]  has  recently  shown  that  the  additional  temperature  dependences  can 
be  understood  in  terms  of  large  variations  in  state-specific  rate  constants,  k,.  In  other  words,  the 
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vibronic  states  of  an  electronic  ground  state  reactant  can  react  at  different  rates.  Birely  and  Lyman  [37] 
have  reviewed  the  state-specific  (or  microscopic)  rate  constant  literature  and  have  found  a  large  body  of 
experimental  evidence  that  points  to  a  trend  towards  greater  reactivity  with  increasing  internal  excita¬ 
tion.  As  an  example.  Spencer  et  al.  [38]  have  considered  the  reaction 

HO  +  H2  —  H20  +  H  (4) 

and  found  that  vibrational  excitation  of  HO  did  not  substantially  affect  the  overall  rate  of  reaction. 
However,  Zellner  et  al.  [28]  observed  that  H2  vibrational  excitation  provided  significant  state-specific 
rate  constant  enhancement.  Specifically,  they  found 

^ho+h,(v- d/^ho+Hj(v-o)  ”  1-5  *  102-  *5) 

Thus,  the  additional  temperature  dependence  of  the  macroscopic  rate  constant  for  the  HO  +  H2  reac¬ 
tion,  as  shown  by  the  concave-upward  curvature  of  the  log  k  ( T)  versus  1000/T  plot  in  Figure  5,  results 
from  the  variation  in  reactivity  of  one  of  its  reactants,  H2(v),  as  a  function  of  temperature.  This  does 
not  require  nonequilibrium  thermal  distributions  of  energy  among  reactants,  but  is  compatible  with  a 
Boltzmann  distribution  of  reactants  over  accessible  vibronic  states. 

Of  exceptional  value  in  the  assembly  of  the  H2-02  reaction  mechanism  has  been  the  compilation 
of  Baulch  et  al.  [39]  at  Leeds  University.  Although  we  conducted  a  literature  search  to  update  the  rate 
constants  in  this  compilation,  very  little  difference  was  found  in  the  rate  constant  reports  following  their 
publication  date  of  1971. 

III.  INDUCTION  TIME  CALCULATIONS 

A.  Background 

Reaction  induction  times,  also  called  ignition  delay  times,  have  become  important  parameters 
measured  in  chemical  kinetic  studies.  Although  a  number  of  workers  had  previously  studied  induction 
times  in  the  H2-02  system  [40-41],  Schott  and  Kinsey  were  one  of  the  first  to  use  these  measurements 


to  examine  H1-O1  reaction  Kinetics  [42].  In  these  experiments,  they  showed  the  inverse  dependence 
of  induction  time  on  [02]  and  were  able  to  calculate  a  rate  coefficient  for  the  important  chain  branching 
reaction,  H  +  02  —  HO  +  O.  Following  these  initial  investigations,  a  large  amount  of  work  was 
devoted  to  the  elucidation  of  the  details  of  the  H2— 02  mechanism  using  induction  time  measurements. 
This  information  has  been  reviewed  by  Schott  and  Getzinger  [43]. 

The  induction  time  of  a  reaction  is  defined  to  be  the  length  of  elapsed  time  between  the  time  of 
reaction  initiation  and  the  time  when  reaction  is  first  observed.  Ambiguity  arises  when  this  definition  is 
applied  to:  1)  a  system  in  which  reaction  initiation  is  not  instantaneous  and  2)  a  system  where  reaction 
observation  depends  upon  the  sensitivity  and  specificity  of  the  observing  apparatus. 

Ambiguity  concerning  reaction  initiation  is  usually  resolved  by  limiting  the  determination  of  reac¬ 
tion  induction  times  to  experimental  apparatus  in  which  the  rise  time  of  a  temperature  jump  or  of  a 
reactive  species  concentration  is  of  the  order  of  a  few  microseconds  or  less.  Reaction  perturbations 
with  rise  times  this  short  have  generally  been  found  to  satisfy  effectively  the  constraint  of  instantaneous 
reaction  initiation. 

The  end  of  the  induction  period  is  often  obvious  since,  even  for  large  variations  in  reaction  condi¬ 
tions,  it  is  usually  associated  with  an  exponential  growth  of  intermediate  species  concentrations.  This 
situation  is  generally  accompanied  by  similar  increases  in  pressure  and  temperature.  Thus,  variations  in 
sensitivity  and  specificity  between  apparatus  conforming  to  the  instantaneous  reaction  initiation  criterion 
should  only  result  in  small  systematic  discrepancies  in  data. 

Because  of  the  above  considerations,  shock-tube  experiments  represent  one  of  the  principal 
methods  used  to  measure  induction  times.  In  addition  to  approximating  instantaneous  reaction  initia¬ 
tion.  these  experiments  usually  provide  at  least  two  other  desirable  features:  1)  isothermal  and  2) 
homogeneous  reaction  conditions. 
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A  near  isothermal  reaction  condition  is  achieved  when  thermochemical  effects  are  mitigated  by 
the  dilution  of  reactive  mixtures.  Diluents  commonly  used  are  inert  gases  such  as  argon  and  helium. 
There  is,  however,  another  factor  peculiar  to  shock-tubes  that  causes  deviations  from  the  ideal  isother¬ 
mal  reaction  condition.  This  has  been  discussed  by  Belles  and  Brabbs  [44]  and  Mirels  [45]  who  have 
shown  that  boundary  layer  growth  in  the  flowing  gases  behind  incident  shocks  can  increase  gas  tempera¬ 
ture,  density  and  particle  residence  time.  In  those  shock-tube  experiments  that  do  not  account  for  this 
effect,  incorrect  temperature  estimates  can  result.  Studies  that  are  most  likely  affected  involve  reactions 
having  high  activation  energies.  Consequently,  it  is  generally  held  that  reaction  rate  coefficients 
obtained  by  shock-tube  measurements  prior  to  —1970  (when  it  became  common  practice  to  make 
boundary  layer  corrections)  should  be  suspected  of  high  apparent  values.  Similarly,  induction  time 
measurements  made  in  shock-tube  experiments  not  correcting  for  boundary  layer  effects  should  be 
suspected  of  being  underestimated. 

Another  advantage  of  shock-tube  experiments  is  that  they  permit  data  to  be  collected  in  time 
periods  that  are  short  in  comparison  to  gas  diffusion  rates.  This  eliminates  the  possibility  of  hetero¬ 
geneous  wall  reactions.  Since  the  elementary  reaction  mechanism  contained  in  Table  I  consists  solely 
of  homogeneous  reactions,  shock-tube  measurements  should  provide  excellent  data  for  testing  numeri¬ 
cal  calculations  using  this  reaction  set. 

Studies  of  the  variation  of  H2-02  reaction  induction  time  as  a  function  of  pressure,  temperature 
and  reactant  composition  have  identified  two  principal  regions  of  ignition.  Schott  and  Kinsey  [42] 
observed  that  above  —  1100°K  and  below  —2  atm  (the  high-temperature,  low-pressure  region,  which 
we  shall  designate  "HTLP")  their  induction  time  measurements  could  be  related  by: 

log|0![O2],r,!  =  A  +  B/T  (6) 

where  [02],  is  the  initial  02  concentration,  r,  is  the  induction  time,  A  and  B  are  constants  and  7"  is  the 
temperature.  They  derived  this  expression  from  a  partial  steady  state  analysis  of  the  following  mechan¬ 
ism: 
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H2.  0;  +  M  —  H.O.HO  +  M 
H  +  02  —  HO  +  0 
0  +  H2  —  HO  +  H 
HO  +  H2  — *  H20  +  H 
H  +  02  +  M  —  H0:  +  M 

Miyama  and  Takeyama  [46]  carried  out  similar  experiments  to  those  of  Schott  and  Kinsey  and 
were  also  able  to  identify  an  inverse  dependence  of  induction  time  on  [02],.  However,  when  Miyama 
and  Takeyama  extended  their  induction  time  experiments  to  temperatures  below  —1100°K  and  pres¬ 
sures  above  ~2  atm  (the  low-temperature,  high-pressure  region,  which  we  shall  designate  "LTHP"). 
they  observed  two  consecutive  temporal  jumps  in  HO  concentration.  The  first  HO  concentration  rise 
was  not  accompanied  by  a  pressure  increase,  whereas  the  second  HO  concentration  was  associated  with 
a  simultaneous  jump  in  pressure.  Their  attempts  to  understand  their  pressure-jump  data  did  not  show  a 
correlation  of  induction  times  with  [02],,  but  did  show  a  correlation  with  [H2],.  Similar  results  were 
obtained  by  Strehlow  and  Cohen  [47],  Fujimoto  [48]  and  Skinner  and  Ringrose  [49],  In  addition  to 
these  observations.  Soloukin  and  coworkers  [50]  and  Meyer  and  Oppenheim  [51]  have  reported  that 
their  schlieren  experiments  have  shown  two  distinct  types  of  ignition.  A  sharp,  or  single  source,  fast 
ignition  was  observed  in  the  HTLP  region  and  a  mild,  or  multiple  source,  slow  ignition  was  observed  in 
the  LTHP  region. 

A  qualitative  understanding  of  the  characteristics  that  differentiate  the  two  regions  of  H2-02  igni¬ 
tion  can  be  obtained  by  considering  the  extended  second  explosion  limit  shown  in  Figure  6.  Brokaw 
[521  has  derived  an  analytic  solution  of  a  simplified  reaction  mechanism  that  simulates  the  different 
behavior  found  in  the  "sharp"  and  "mild"  ignition  region.  In  addition,  Wakefield,  Ripley  and  Gardiner 
[53]  have  demonstrated  the  origins  of  some  of  those  differences  by  comparing  analytic  and  numeric 
solutions  of  a  more  detailed  reaction  mechanism.  In  brief,  the  results  of  these  studies  have  shown  that 
the  two  regions  of  ignition  are  separated  by  the  extended  second  explosion  limit,  at  which  the  rate  of 
chain  branching  is  balanced  by  the  rate  of  chain  termination.  To  the  left  of  the  extended  second  limit, 
in  the  LTHP  region,  the  H02  concentration  becomes  significant  and  its  bimolecular  reactions  become 
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important  in  mild  ignition.  To  the  right  of  the  extended  second  limit,  in  the  HTLP  region,  chain 
branching  reactions  among  H,  0  and  HO  dominate  and  lead  to  sharp  ignition. 

B.  Results 

The  information  discussed  above  indicates  that  homogeneous  reactions  predominate  in  both  the 
LTHP  and  HTLP  regions  of  ignition  in  H2-02  gas  mixtures.  Thus,  we  have  calculated  H2-02  reac¬ 
tion  induction  times  for  both  regions  using  the  mechanism  given  in  Table  I.  For  the  LTHP  region,  we 
have  computed  a  set  of  induction  times  for  initial  conditions  descriptive  of  the  work  by  Skinner  and 
Ringrose  [49],  These  are  shown  in  Figure  7.  For  the  HTLP  region,  we  have  calculated  a  series  of 
induction  times  using  the  reaction  conditions  of  White  [54] .  These  data  are  presented  in  Table  3.  The 
agreement  between  the  calculated  and  experimental  induction  times  found  in  both  the  HTLP  and 
LTHP  regions  of  H2-02  ignition  supports  the  mechanism  of  Table  1  in  its  description  of  the  initial 
phases  of  reaction. 

As  a  means  of  investigating  the  functional  behavior  of  induction  time  in  the  H2-02  system,  we 
have  performed  a  number  of  adiabatic  and  isothermal  calculations  while  varying  initial  temperature, 
pressure  and  extent  of  ,V2  dilution.  Under  the  conditions  of  these  calculations,  <V2  does  not  participate 
as  a  reactant  but  only  as  an  inert  diluent.  We  show  the  adiabatic  temperature  versus  reaction  time 
curves  that  were  generated  for  a  stoichiometric  H2-air  mixture  in  Figures  8-11.  In  these  curves,  the  ini¬ 
tial  temperature  was  varied  from  900  to  2000°K  in  100°K  increments  for  each  of  the  initial  pressures 
0.05,  0.5,  5  and  50  atm.  The  times  at  which  the  H  atom  and  HO  radical  intermediates  experience  their 
concentration  maxima  appear  to  coincide  and  are  indicated  by  arrows  on  each  curve.  It  is  noted  that 
these  times  also  correspond  to  the  first  temperature  increases  and  exponentially  increase  with  decreasing 
initial  temperature.  It  is  also  seen  that  the  curves  generated  in  the  HTLP  regions  show  a  distinctly 
different  form  than  the  curves  generated  in  the  LTHP  regions.  These  observations  seem  to  agree  with 
the  different  ignition  patterns  reported  for  these  two  regions. 
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In  Figures  12-14.  we  show  the  induction  time  dependence  on  initial  temperature  and  pressure  as  a 
function  of  dilution  of  a  stoichiometric  Hi  —  Oi  mixture.  The  change  in  ignition  pattern,  as  demon¬ 
strated  by  induction  time,  between  the  LTHP  and  HTLP  regions  is  observed  as  a  shift  toward  sharper, 
exponential  increases  in  induction  time  with  decreasing  initial  temperature  and  increasing  initial  pres¬ 
sure.  The  dilution  effect  appears  to  uniformly  increase  induction  time  for  the  stoichiometric,  reactive 
mixtures  and  does  not  significantly  change  the  form  of  the  induction  time  dependence  on  initial  tem¬ 
perature  and  pressure.  These  observations  are  in  agreement  with  those  of  Meyer  and  Oppenheim  [51]. 

The  change  in  behavior  of  induction  time  and,  thus,  ignition  may  possibly  be  better  demonstrated 
in  Figure  15.  In  this  figure,  induction  time  is  plotted  as  a  function  of  initial  pressure  for  each  of  three 
isothermal  calculations.  Initial  temperature  conditions  were  1000,  1100,  and  1500°K.  It  can  be  seen 
that  the  pressure  dependence  of  induction  time  clearly  makes  a  shift  below  1100°K  and  above  —2  atm. 
Mild  ignition  affected  by  hydrodynamic  interactions  with  chemistry  can  be  easily  imagined  in  a  situation 
where  induction  limes  are  relatively  long. 

IV.  SECOND  EXPLOSION  LIMIT  CALCULATIONS 

A.  Background 

The  line  showing  the  division  between  explosive  and  nonexplosive  regions  (fast  and  slow  reac¬ 
tion)  for  a  stoichiometric  Hi-02  mixture  contained  in  a  KCl-coated  reaction  vessel  has  been  shown  in 
Figure  1.  The  line  makes  its  excursion  between  — 650-850°K  and  between  — 1-5000  torr  and  is  divided 
into  three  segments  designated  as  the  first,  second  and  third  explosion  limits. 

As  mentioned  previously,  the  three  explosion  limits  have  been  investigated  in  static  flow  experi¬ 
ments  that  have  demonstrated  how  important  the  dimensions  and  materials  of  the  reaction  vessels  are 
in  determining  points  of  explosion.  These  observations  firmly  established  the  significance  of  gas- 
dynamic  processes  in  these  experiments.  They  have  also  shown  that  the  first  and  second  explosion  lim¬ 
its  depend  upon  the  chain  branching  chemistry  of  the  H2-Oj  system  [22].  However,  there  is  still 
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some  question  concerning  the  mechanism  of  the  third  explosion  limit,  i  e..  whether  it  is  controlled  by 
thermal  or  chain  branching  processes  [551 . 


In  order  to  describe  the  influence  of  chain  reactions  at  the  first  and  second  explosion  limits,  it  has 
become  common  practice  to  consider  the  concentration  of  a  typical  chain  carrier  [431 

J[C\IJt  -  /  +  fl(C|  -  HC1  -  /  +  V(C'|  (7) 

where  [Cl  represents  chain  carrier  concentration.  /  is  the  initiation  rate.  B  is  the  rate  of  chain  branching 
and  7"  is  the  chain  termination  rate.  The  V  of  Equation  (7)  is  the  net  sum  of  the  rates  of  chain  branch¬ 
ing  and  chain  termination  and  is  an  important  parameter  used  in  describing  the  concentration  of  a  chain 
carrier  in  the  early  stages  of  reaction. 

Assuming  initial  reaction  rate  conditions,  where  /.  B  and  T  are  well-delined  and  [C]  is  negligible. 
Equation  (7)  may  be  integrated  to  give 

[Cl  -  (//.V)  (exp(.W)  -  I|.  (8) 

Two  limits  of  this  equation  describe  the  chemistry  on  opposite  sides  (slow  versus  fast  reaction)  of  an 
explosion  limit  controlled  by  chain  reaction  chemical  kinetics  In  the  first  case,  where  chain  termina¬ 
tion  dominates  chain  branching  li  e..  T  >  B),  \  is  negative  and  [Cl  is  forced  to  approach  a  concentra¬ 
tion  approximated  by  <//.V)  In  the  second  case,  where  chain  branching  overshadows  chain  termina¬ 
tion  <i.e..  B  >  D.  .V  is  positive  and  [Cl  exponentially  increases  with  time.  In  this  latter  situation,  a 
large  amount  of  heat  is  released  in  a  short  time  interval  and  the  total  complex  of  elementary  reactions 
is  accelerated  to  fast  reaction  (explosion)  by  thermal  feedback  or  autocatalysts. 


The  Hj-0:  mechanism  responsible  for  this  behavior  may  be  summarized  in  terms  of  its  most 
important  reactions.  The  three  elementary  reactions  of  major  significance  to  chain  branching  in  the 
vicinity  of  the  first  and  second  explosion  limns  are  [421: 

H+O;— HO  +  O  (9 > 

O  H:  —  HO  *  H 
HO  +  H;  -  H  +  H:0 


1 


13 


(10) 

(ID 
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Similarly,  chain  mutation  involves  the  four-center  reaction  [56]: 

H;  +  0, -2HO  (12) 

although  a  number  of  alternate  initiation  reactions  may  play  an  important  role.  Where  the  reactivity  of 
H  >  >  HO-,  chain  termination  occurs  principally  by  the  reactton  [13]: 

H  +  0:  +  U  -  H0:  +  Xi.  (13) 

In  addition,  two  other  termolecular  reactions  have  been  shown  to  be  important  chain  termination  reac¬ 
tions  [57] 

H  +  HO  +  W  -  H:0  +  M  (14) 

H  +  H  +  .V#  —  H,  +  Xi  (15) 

Although  these  elementary  reactions  do  not  constitute  a  complete  reaction  mechanism,  they  are  pri¬ 
marily  responsible  for  many  of  the  features  of  the  first  and  second  explosion  limits. 

B.  Results 

Because  the  second  explosion  limn  is  largely  determined  by  the  homogeneous  chemical  kinetics  of 
the  H;-0-  system,  we  have  tested  our  mechanism  by  performing  a  number  of  calculations  of 
parameters  associated  with  this  limit  Although  several  other  numerical  studies,  using  either  abbrevi¬ 
ated  reaction  mechanisms  or  varying  individual  reaction  rate  coefficients,  have  examined  the  second 
explosion  limit  [55.58.59],  we  have  based  our  evaluations  on  straightforward  computations  involving  a 
complete  reaction  set  We  have  not  varied  the  rate  coefficients  to  achieve  better  fits  to  data.  In  these 
calculations,  it  was  not  our  purpose  to  fit  a  parameter  but  rather  to  compare  computed  and  experimental 
results  and  to  try  to  understand  the  significance  of  any  agreement  or  disagreement. 

As  a  first  step,  we  have  made  a  series  of  calculations  describing  a  0  5  atm.  H2:0;/2: 1  reaction  mix¬ 
ture  at  10°K  intervals  from  808  to  848°K  The  behavior  of  intermediate  species  concentrations  was  fol¬ 
lowed  temporally  and  was  recorded  as  a  function  of  initial  temperature 
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As  shown  in  Figures  12-14,  we  found  that  reaction  induction  time  increases  exponentially  with 
decreasing  initial  temperature.  This  behavior  is  also  observed  in  the  temporal  plots  of  the  hydroxyl  rad¬ 
ical  concentration  presented  in  Figure  16.  The  pattern  continues  up  to  a  temperature  between  818  and 
828°K.  where  further  reduction  of  initial  temperature  results  only  in  linear  increases  in  induction  time. 
At  about  the  same  temperature,  a  part  of  the  exponential  growth  in  HO  concentration  appears  to  van¬ 
ish.  This  transition  occurs  within  10°K  and  appears  to  be  linked  to  the  loss  of  ignition  at  the  second 
explosion  limit.  In  addition,  we  believe  that  this  latter  feature  may  also  be  tied  to  the  experimental 
observation  of  spikes  and  overshoots  of  intermediate  species  concentrations  noted  in  shock-tube  experi¬ 
ments  conducted  close  to  the  second  explosion  limit  [43], 

Using  the  above  changes  in  induction  time  and  hydroxyl  radical  behavior  to  define  the  crossing  of 
the  second  explosion  limit,  we  have  computed  a  set  of  points  in  the  temperature-pressure  plane  that 
corresponds  to  the  second  explosion  limit  of  a  stoichiometric  H:-0;  gas  mixture.  These  results  are 
presented  in  Figure  17  where  they  are  compared  to  the  experimental  work  of  Lewis  and  von  Elbe  [11]. 
We  note  that  the  calculated  values  show  a  second  explosion  limit  uniformly  lower  in  temperature  than 
that  which  was  observed  in  a  7  4  cm  diameter.  KCI-coated  reaction  bulb.  This  discrepancy  can  be  attri¬ 
buted  to  the  fact  that  we  have  modeled  homogeneous  chemistry  in  a  situation  where  heterogeneous 
chemistry  also  plays  a  part.  The  displacement  of  our  calculated  data  towards  lower  temperatures  is 
consistent  with  the  idea  that  the  absence  of  heterogeneous  termination  reactions  in  our  mechanism 
would  require  less  thermal  energy  supplied  to  branching  reactions  to  overbalance  chain  termination. 

By  examining  the  production  and  loss  terms  for  each  chemical  species,  we  have  observed  that  the 
major  reaction  channel  for  consumption  of  HO  changes  as  the  initial  reaction  temperature  is  reduced 
below  the  second  explosion  limit.  Specifically,  the  chain  branching  reaction.  H  +  0:  — •  HO  +  0 
(Equation  9),  is  replaced  by  the  chain  termination  reaction.  H  4-  0:  +  M  —  HO:  +  W  (Equation  13) 
This  observation  appears  to  coincide  with  the  disappearance  of  the  second  [HO]  exponential  increase 
shown  in  Figure  16  On  the  explosive  side  of  the  second  explosion  limit  and  at  the  point  of  ignition. 
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reaction  (13)  becomes  an  insignificant  channel  for  production  of  H02  and  reaction  (9)  becomes  a  major 
channel  for  production  of  HO.  On  the  nonexplosive  side  of  the  second  explosion  limit,  reaction  (13) 
remains  the  primary  channel  of  formation  of  H02,  with  an  increase  in  H02  concentration  followed  by 
enhanced  reaction  rates  for  H02  and  H202. 

Defining  the  initial  reaction  rate  to  be  proportional  to  the  rate  of  formation  of  H20,  we  have  com¬ 
puted  this  parameter  above  the  second  explosion  limit.  Figure  18  presents  the  results  of  a  calculaltion 
of  </[H20]/dt  as  a  function  of  pressure  along  a  760°K  isotherm  for  a  stoichiometric  H2-02  gas  mixture. 
The  behavior  of  the  calculated  initial  reaction  rate  above  the  second  explosion  limit  closely  resembles 
that  experimentally  determined  by  Lewis  and  von  Elbe  [60].  As  the  pressure  of  the  reactive  gas  mix¬ 
ture  is  raised  above  the  second  explosion  limit,  the  initial  reaction  rate  falls  rapidly  from  extremely 
large  values  to  a  minimum  and  then  begins  to  gradually  climb  as  pressure  is  increased. 

V.  THE  TEMPORAL  BEHAVIOR  OF  REACTION  SPECIES 

A.  Background 

Since  a  good  deal  is  known  about  the  individual  elementary  reactions  of  the  H2-02  system 
(specifically,  their  rate  coefficients  as  a  function  of  temperature),  it  is  possible  to  test  the  complete  reac¬ 
tion  mechanism  by  comparing  computed  and  observed  temporal  profiles  of  reaction  species  concentra¬ 
tions.  In  the  shock-tube,  temporal  profiles  of  species  concentrations  have  been  determined  typically  as 
a  function  of  the  stage  of  H2  combustion.  Moreover,  these  experimental  profiles  have  not  been  pub¬ 
lished  in  simple,  time-based  plots  of  species  concentrations,  but  are  commonly  used  to  derive  other 
quantities  such  as  exponential  growth  parameters,  which  are  discussed  below. 

The  H2-02  reaction  system,  under  low-density  conditions,  has  been  shown  to  proceed  by  well- 
defined  stages  of  combustion:  i.e..  1)  induction.  2)  transition,  3)  partial  equilibrium,  and  4)  full  equili¬ 
brium.  Chain  initiation,  branching  and  propagation  reactions  appear  to  dominate  the  induction  stage  of 
combustion.  The  transition  stage  primarily  involves  chain  branching  and  propagation  and  is  the  region 
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wherein  species  concentration  spikes  and  overshoots  are  observed.  Association  reactions  become  impor¬ 
tant  in  the  partial  equilibrium  stage  of  system  reaction,  and.  finally,  the  fuil  equilibrium  stage  of 
combustion  occurs  when  all  reaction  rates  become  balanced. 

B.  Results 

Since  the  hydroxyl  radical  is  the  intermediate  species  most  often  studied  in  shock-tube  investiga¬ 
tions  of  the  Hi  — O;  system  [43|,  we  have  calculated  its  behavior  in  the  various  stages  of  Hi  combus¬ 
tion.  Generally  speaking,  direct  observation  of  the  other  intermediate  species  has  been  hampered  by  a 
lack  of  methods  having  the  required  sensitivities  and  response  times.  The  O  atom  concentrations  have 
been  measured,  but  primarily  by  methods  that  involve  the  presence  of  CO  as  a  reactant  [61],  This 
technique  requires  another  set  of  elementary  reactions  that  must  be  considered  before  the  O  atom  con¬ 
centrations  can  be  calculated.  The  H  atom  concentrations  have  recently  been  measured,  but  these 
results  appear  to  suffer  from  the  presence  of  impurities  (62].  The  temporal  behavior  of  HOi  in  shock- 
tube  work  is  not  well  known. 

Throughout  most  of  the  induction  stage,  the  HO  concentration  has  an  exponential  growth  rate  and 
may  be  expressed  as: 

[HO]  —  (HOI,,  exp(«r )  (16) 

where  r  is  time  and  <  is  the  exponential  growth  parameter  that  depends  upon  chain  branching,  propaga¬ 
tion  and  termination  as  outlined  in  Section  IV.  Jachimowski  and  Houghton  (63)  have  investigated 
[HOI  in  the  induction  period  using  incident  shock-wave  experiments  coupled  with  ultraviolet  absorption 
spectroscopic  detection  of  HO.  From  their  data,  they  calculated  «  for  a  variety  of  fuel-rich  and  fuel- 
lean  HiOi  Ar  gas  mixture  at  different  temperatures.  Typically,  they  identified  an  induction  time 
corresponding  with  a  detectable  level  of  HO  concentration  and  defined  [HOI,,  as  a  pseudo-initial  con¬ 
centration  related  to  the  rate  of  initiation.  By  rewriting  Equation  (16)  as. 

r«  -  In  [HOI./IHOI..  (17) 
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they  were  able  to  calculate  the  pseudo-initial  hydroxyl  concentration.  [HO]„,  as  a  function  of  their 
graphically  measured  induction  time,  r,  exponential  growth  factor,  e,  and  induction-time-related  HO 
concentration,  [HO].. 

We  have  performed  several  computations  involving  the  induction  stage  of  combustion  using 
Jachimowski  and  Houghton’s  initial  reaction  conditions.  In  Table  4.  we  present  a  comparison  of  their 
experimental  and  our  calculated,  [HO]-based,  exponential  growth  parameters.  A  comparison  of  «/[0;]„ 
versus  1000/T  is  also  provided  in  Figure  19  to  demonstrate  the  temperature  and  reactant  concentration 
dependencies  of  the  HO  exponential  growth  parameter. 

Getzinger  and  Schott  have  examined  the  partial  equilibrium  stage  of  H>  combustion  using  HO 
concentration  profiles  [64],  In  their  experiments,  they  initiated  the  reaction  using  incident  shock-waves 
and  then  followed  the  progress  of  the  reaction  using  HO  ultraviolet  absorption  spectroscopy.  The  par¬ 
tial  equilibrium  stage  of  combustion  occurs  when  the  reaction  cycle  reaches  a  point  at  which  the  ter- 
molecular  association  reactions  decelerate  the  overall  reaction  rate.  The  result  is  that  the  reaction 
species  concentrations  remain  essentially  constant  for  a  protracted  time  interval.  These  concentrations 
may  be  calculated  from  a  single  known  concentration  once  the  appropriate  equilibrium  constant  infor¬ 
mation  is  known. 

Following  the  experimental  measurement  of  the  [HO]  profiles  for  a  number  of  fuel-lean  H:-0: 
mixtures,  Getzinger  and  Schott  calculated  the  temporal  behavior  of  the  remaining  reaction  species  using 
analytic  functions  describing  the  partial  equilibrium  condition.  In  order  to  focus  on  the  recombination 
zone  of  reaction,  they  calculated  concentrations  in  terms  of  a  dimensionless  reaction  progress  variable, 
y,  defined  by  the  equation: 

y  -  (N  -  Sfq)H\  -  \eq)  (18) 

where  ,V  is  the  instantaneous  species  concentration  and  .V,,,  is  the  species  concentration  at  full- 
equilibrium.  Thus,  reaction  initiation  occurs  at  y  “  1  and  full-equilibrium  is  attained  at  y  —  0. 
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We  have  calculated  the  temporal  behavior  of  each  of  the  reaction  species  concentrations  as  a  func¬ 
tion  of  y  for  one  of  the  fuel-lean  Hi  —  02  mixtures  examined  by  Getzinger  and  Schott.  Our  data  is 
presented  in  Figure  20  and  includes  a  set  of  points  representative  of  those  obtained  experimentally  for 
[HOI.  The  calculated  data  appears  to  be  almost  identical  with  the  measured  data. 

In  addition,  we  have  calculated  the  temporal  behavior  of  the  reaction  species  number  densities  for 
the  same  Hi.Oi.Ar  mixture  used  above  and  for  a  mixture  containing  no  Ar.  We  show  these  data  in 
Figures  21  and  22.  These  figures  demonstrate  the  increased  reaction  time  effect  that  dilution  has  on 
the  partial  equilibrium  stage  of  combustion.  They  also  show  the  relation  of  the  partial  equilibrium  stage 
of  reaction  to  the  other  stages  of  combustion. 

VI.  REACTION  EQUILIBRIUM  CALCULATIONS 

Another  test  of  the  detailed  H2-0:  mechanism  is  to  determine  whether  the  mechanism  can  predict 
the  species  concentrations  once  chemical  equilibrium  is  reached.  Under  adiabatic  conditions,  this 
means  that  certain  stale  and  thermodynamic  variables  are  accurately  determined.  In  this  section,  we 
compare  the  long-time  equilibrium  results  of  our  computations  with  those  obtained  using  the  NASA- 
LEWIS  "Computer  Program  for  Calculation  of  Complex  Chemical  Equilibrium  Compositions"  [65]. 

The  comparison  described  in  this  section  is  between  two  entirely  different  computational  methods. 
The  NASA-LEWIS  Program  uses  a  minimization  of  Gibbs  (or  Helmholtz)  free-energy  technique  to 
derive  the  composition  of  the  chemical  equilibrium  mixture,  whereas  our  calculations  involve  the 
numerical  integration  of  time-dependent,  ordinary  differential  equations  defined  by  the  reaction 
mechanism  contained  within  Table  1.  The  use  of  a  search  routine  by  the  NASA-LEWIS  Program  to 
define  the  reaction  species  included  in  their  calculation  adds  generality  to  the  comparison. 

We  have  carried  out  calculations  with  both  programs  to  obtain  the  chemical  equilibrium  composi¬ 
tion  of  an  adiabatic,  constant-volume  combustion  of  one  atmosphere  of  stoichiometric  hydrogen-air 
< H2:0::N;/ 2:1:4)  with  an  initial  temperature  of  900°K.  A  comparison  of  our  computed  data  is 
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presented  in  Table  5.  In  addition  to  the  reaction  species  concentrations,  the  chemical  equilibrium  mix¬ 
ture  density,  pressure  and  temperature  are  calculated.  The  excellent  correlation  between  the  two  sets  of 
computed  data  substantiates  the  detailed  mechanism  of  Table  1  as  a  good  description  of  the  reactive 
system  in  its  approach  to  chemical  equilibrium. 

Moreover,  the  NASA-LEWIS  Program  calculations  are  in  agreement  with  our  initial  assumptions 
regarding:  (1)  the  inert  character  of  N2asa  diluent,  and  (2)  the  predominant  role  of  neutral  species  in 
the  H2  —  02  reactive  system.  Several  NASA-LEWIS  Program  computations  were  made  in  which  ionic 
species  and  nitrogen-containing  reaction  species  were  allowed  to  participate  in  the  chemical  equilibrium 
mixture.  In  all  situations  where  the  final  adiabatic  temperature  was  below  ~3000°K,  it  was  found  that 
;V2  remained  essentially  inert  and  that  ionic  species  concentrations  were  negligible. 

VII.  SUMMARY  AND  CONCLUSIONS 

The  principal  goal  of  the  work  presented  in  this  paper  was  to  determine  the  accuracy  with  which  a 
proposed  reaction  mechanism  is  able  to  model  the  chemistry  of  the  H2-0:  reactive  system.  Since 
analytical  and  intuitive  interpretations  of  experimental  evidence,  coupled  with  conjecture  and  supposi¬ 
tion,  are  the  common  tools  used  in  the  construction  of  a  detailed  reaction  mechanism,  we  have 
searched  for  experimental  data  relatively  free  of  hydrodynamic  and  diffuse  transport  effects  against 
which  we  could  test  the  proposed  mechanism.  Although  it  is  not  possible  to  define  measurable  parame¬ 
ters  that  are  unaffected  by  gasdynamic  interactions,  some  that  are  minimally  affected  include:  (1)  the 
reaction  induction  (delay,  ignition)  times,  (2)  the  temperature-pressure  dependence  of  the  second 
explosion  limit,  (3)  the  rate  of  reaction  above  the  second  explosion  limit,  (4)  the  temporal  behavior  of 
reaction  species,  under  certain  circumstances  (i.e.,  shock  tubes). 

The  elementary  reactions  most  important  in  determining  the  ignition  properties  of  the  H2-02 
system  are  those  involved  in  chain  initiation,  branching  and  propagation.  Reaction  induction  times, 
which  depend  on  the  integrated  effect  of  these  reactions,  have  served  as  a  good  test  of  the  mechanism 
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in  this  regime.  The  agreement  between  experimental  and  computational  ignition  delays  appears  to 
confirm  the  set  of  elementary  reactions  in  its  description  of  the  initial  reaction  phase  of  the  H2  — 02  sys¬ 
tem. 


More  specifically,  we  have  shown  that  calculations  in  the  high-temperature,  low-pressure  (HTLP) 
region  of  H2  combustion  correlate  well  with  experimental  induction  times.  This  region  is  characterized 
by  strong,  single-source  or  integrated  kernel  ignition  having  a  sharply  defined  blast-front.  The  sharp¬ 
ness  and  uniformity  of  this  strong  ignition  implies  that  timescales  for  chemical  interactions  are  much 
shorter  than  timescales  for  hydrodynamic  interactions  and,  thus,  we  might  expect  good  agreement 
between  the  calculations  and  experiment. 

However,  ignition  in  the  low-temperature,  high-pressure  (LTHP)  region  is  affected  by  the  cou¬ 
pling  of  gasdynamics  and  chemistry.  This  is  demonstrated  by  the  mild,  multiple  source  or  multiple- 
kernel  appearance  of  laser-schlieren  experiments.  Even  though  this  region  represents  the  worst  case  for 
zero-dimensional  calculations  of  induction  times,  there  is  good  agreement  between  computed  and 
experimental  induction  times. 

The  second  explosion  limit  calculations  provide  a  different  kind  of  test  of  the  H2  — 02  mechanism. 
In  these  calculations,  third-body  recombination  or  association  reactions  are  important  and  compete  with 
the  bimolecular  chain-branching  and  propagation  reactions.  Our  calculations  were  performed  for 
stoichiometric  H2-02  mixtures  in  the  slow-reactive  zones  of  the  temperature-pressure  plane.  They 
show  that,  for  a  given  pressure,  the  calculated  and  measured  explosion  limits  disagree  by  about  20°K. 
This  discrepancy,  however,  is  understandable  when  we  take  into  account  the  fact  that  the  first  and 
second  explosion  limits  are  affected  by  vessel  dimensions  and  materials.  In  the  experimental  results 
shown  in  Figure  17,  there  is  a  contribution  to  the  exponential  growth  factor,  «,  of  Equation  6  due  to 
chain  termination  by  heterogeneous  reactions  on  the  chamber  walls.  Thus,  the  need  to  include  hetero¬ 
geneous  reactions  indicates  an  important  difference  between  our  mechanism  and  one  necessary  for  the 
description  of  this  specific  experiment.  Nevertheless,  these  explosion  limit  calculations  have  shown  the 


H'HWrr-  r»»> 


I 


BURKS  AND  ORAN 

correct  trend  of  the  second  explosion  limit.  The  displacement  of  the  calculated  results  toward  lower 
temperatures  is  reasonable  when  the  absence  of  heterogeneous  reactions  in  our  mechanism  is  con¬ 
sidered. 

An  additional  test  of  the  mechanism  in  Table  1  was  provided  by  computations  carried  out  along 
an  isotherm  in  the  slow-reactive  zone  above  the  second  explosion  limit.  The  appropriate  balance 
between  bimolecular  chain  branching  and  termolecular  chain  termination  is  demonstrated  by  the  calcu¬ 
lated  results  for  the  overall  reaction  rate  producing  H20  as  presented  in  Figure  18. 

Comparisons  were  also  made  between  calculations  and  experimental  measurements  of  the  tem¬ 
poral  behavior  of  hydroxyl  radicals  under  different  reaction  conditions  found  in  shock-tube  experi¬ 
ments.  These  HO  concentrations  were  first  calculated  for  the  induction  period  of  highly-diluted. 
Hi  —  02  gas  mixtures.  Exponential  growth  factors,  e.  were  calculated  for  a  number  of  temperatures, 
pressures,  and  Hi-Oi-dilueni  mixture  compositions.  As  shown  in  Table  4,  the  experiments  and  calcu¬ 
lations  are  in  good  agreement.  These  tests  also  help  substantiate  the  chain  initiation,  branching  and 
propagation  parts  of  the  reaction  mechanism. 

The  partial  equilibrium  stage  of  Hi  combustion  is  the  region  in  which  termolecular  association 
reactions  become  important  to  the  deceleration  of  chain-linked  reactions.  Here  is  where  the  major 
amount  of  thermochemical  energy  release  is  achieved.  The  temporal  behavior  of  HO  concentrations  in 
this  stage  was  calculated  and  is  shown  together  with  the  experimental  results  in  Figure  20.  Here,  again, 
good  agreement  between  experiment  and  computation  was  obtained.  This,  as  in  the  second  explosion 
limit  calculations,  adds  support  to  the  validity  of  the  chain  termination  reactions. 

Finally,  the  proposed  reaction  mechanism  was  used  to  calculate  the  species  concentrations,  den¬ 
sity,  temperature  and  pressure  of  a  reaction  mixture  that  had  reached  chemical  equilibrium.  The  good 
agreement  obtained  between  the  results  of  this  calculation  and  a  calculation  and  a  calculation  made 
using  the  NASA-LENIG  Program  (which  is  a  technique  for  calculating  equilibrium  properties  based  in 
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minimization  the  the  free  energy)  supports  the  general  validity  of  the  proposed  reaction  mechanism, 
the  reaction  rate  coefficients,  and  the  species  enthalpies. 
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Table  1.  H;  -  O;  Elementary  Reaction  Mechanism. 


— 1 

— 

Reaction 

k.  —  AT8  exp  (- 

-C/T)'J> 

A'”' 

B 

C,b' 

1 

H  +  HO  =  0  +  H; 

1.401-14) 

1.00 

3.501+03) 

Baulch  et  al.,  1972 

-'l 

3.001-14) 

1.00 

4.481+03) 

Baulch  et  al..  1972 

2 

H  +  HO:  =  H;  +  0: 

4.201-11) 

0.00 

3.501  +  02) 

Baulch  et  al.,  1972 

-2 

9.101-11) 

0.00 

2.911+04) 

Baulch  et  al.,  1972 

3 

H  +  HO;  =  HO  +  HO 

4.201-10) 

0.00 

9.501+-2) 

Baulch  et  al.,  1972 

-3 

2.001-11) 

0.00 

2.021  +  04) 

Baulch  et  al.,  1972 

4 

H  +  H;  =  O  +  H;0 

8.301-11) 

0.00 

5.001+02) 

Hampson,  1973 

-4 

1.751-12) 

0.45 

2.841+04) 

kr  -  kf/K, 

5 

H  +  H;0;  =  HO;  f  H; 

2.801-12) 

0.00 

1.901+03) 

Baulch  et  al.,  1972 

~5 

1.201-12) 

0.00 

9.401+03) 

Baulch  et  al.,  1972 

6 

H  +  H;0;  =  HO  +  H;0 

5.281-10) 

0.00 

4.501+03) 

Baulch  et  al.,  1972 

-6 

3.991-10) 

0.00 

4.051  +  04) 

kr  -  kr/K* 

7 

H  +  0,  =  HO  +  0; 

2.701-11) 

0.00 

O.OO(-fOO) 

Hampson,  1973 

-7 

5.201-15) 

0.69 

3.861  +  04) 

k,  -  kf/lCc 

8 

HO  +  H,  =  H  4-  H,0 

1.831-15) 

1.30 

1.841  +  03) 

Cohen  and  Westberg,  1978 

-8 

1.791-14) 

1.20 

9.611+03) 

Cohen  and  Westberg,  1978 

9 

HO  +  HO  =  H;  +  0; 

1.091-13) 

0.26 

1.471  +  04) 

kr  -  k.Kc 

-9 

2.821-11) 

0.00 

2.421  +  04) 

Olson,  1977 

10 

HO  +  HO  rO  +  H;0 

1.001-16) 

1.30 

0.001+00) 

Cohen  and  Westberg,  1978 

-10 

3.201-15) 

1.16 

8.771+03) 

k,  -  kf/Kc 

11 

HO  +  HO;  =  H;0  +  0; 

8.301-11) 

0.00 

5.031  +  02) 

Lloyd,  1974 

-11 

2.381-10) 

0.17 

3.691+04) 

k,  -  kf/Kc 

12 

HO  +  H;0;  =  HO;  +  H; 

1.701-11) 

0.00 

9.101+02) 

Baulch  et  al.,  1972 

-12 

4.701-11) 

0.00 

1.651+04) 

Baulch  et  al.,  1972 

13 

HO  +  0,  =  HO;  4-  0; 

1.601-12) 

0.00 

9.561+02) 

Hampson,  1973 

-13 

6.691-14) 

0.33 

2.041  +  04) 

kr  -  kr/Kc 

14 

HO;  +  H;  =  HO  4-  H;0 

1.201-12) 

0.00 

9.411+03) 

Olson,  1977 

-14 

1.331-14) 

0.43 

3.621  +  04) 

kr  -  kf/Kc 

15 

HO;  4-  HO;  =:  H;0;  4-  0; 

3.001-11) 

0.00 

5.001+02) 

Hampson.  1973 

-15 

1.571-09) 

-0.38 

2.201  +  04) 

k,  -  kf/Kc 
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Table  1  (Continued).  H;-  O?  Elementary  Reaction  Mechanism. 

1  k.  -  AT8  exp  <— C.  T)'J‘  1 


1 

iveawuuii 

— 

A"” 

B 

£»h ) 

16 

HO;  +  Oj  =  HO  +  20; 

1  701-13) 

0  00 

1  411+03) 

Lloyd.  1974 

-16 

9  941-46) 

0.16 

1.60(4-04) 

kr  -  kf/K* 

17 

O  +  HO  =  H  +  02 

2.721-12) 

0.28 

-8  101+01) 

k.  -  k,K. 

-17 

3.701-10) 

0.00 

8  451+03) 

Baulch  et  al..  1972 

18 

O  +  HO;  =2  HO  +  O; 

8.321-11) 

0.00 

5  031+02) 

Lloyd.  1974 

-18 

K* 

K* 

O 

T 

0.18 

2.821+04) 

k.  -  k./K, 

19 

0  +  H;0;  =  H;0  +  0, 

1  401-12) 

0.00 

2.121+03) 

Hampson.  1973 

-19 

5.701-14) 

0.52 

4.481+04) 

k,  -  k,/fC 

20 

0  -  H;0;  =  HO  +  HO; 

1.401-12) 

0.00 

2.131+03) 

Hampson.  1973 

-20 

2.071-15) 

0.64 

8.231  +03) 

k,  -  k./K, 

2! 

0  4-  0?  —  Oj  4*  O2 

1  891-11) 

0.00 

2  301+03) 

Hampson.  1973 

-21 

1.981-11) 

0.00 

5.061  +  04) 

k,  -  kr/K. 

22 

h  +  h  +  m  =  h2  +  m 

1 .80< — 30) 

-1.00 

0.001+00) 

Baulch  et  al..  1972 

-22 

3.701-10) 

0.00 

4.831-04) 

Baulch  et  al..  1972 

23 

H  +  HO  +  M  =  HjO  +  M 

6.201-26) 

-2.00 

0  001+00) 

Baulch  et  al..  1972 

-23 

5.801-09) 

0.00 

5.291+04) 

Baulch  et  al..  1972 

24 

H  +  0;  +  M  r  HOi  +  M 

4.141-33) 

0.00 

-5.001+02) 

Baulch  et  al..  1972 

-24 

3.501-09) 

0.00 

2.301+04) 

Baulch  et  al.,  1972 

25 

HO  +  HO  +  M  =r  HjO;  +  M 

2.501-33) 

0.00 

-2.551  +  03) 

Baulch  et  al..  1972 

*25 

2.001-07) 

0.00 

2.291+04) 

Baulch  et  al..  1972 

26 

0  +  H  +  M  =  H0  +  M 

8.281-29) 

-1.00 

0  001+00) 

Bahn.  1969 

-26 

2.331-10) 

0.21 

5.101+04) 

k,  -  k,/  K. 

27 

0  +  HO  +  M  =  H02  +  M 

2.801-31) 

0.00 

0.001+00) 

Bahn.  1969 

-27 

1.101-04) 

-0.43 

3.221+04) 

kr  -  kr/  K, 

28 

O  +  O  +  MrOi  +  M 

5.201-35) 

0.00 

-9.001+02) 

Baulch  et  al.,  1972 

-28 

3.001-06) 

-1.00 

5.94(4-04) 

Baulch  et  al..  1972 

29 

0  +  02  +  M  =  0,  +  M 

3.311-35) 

0.00 

-9.801  +  02) 

kr  -  k,K,. 

-29 

7.161-10) 

0.00 

1.121+04) 

Heimerl.  1979 

'•■’Bimolecular  reaction  rate  constants  are  given  in  units  of  cm1/ (molecule  sect. 

,bl Exponentials  to  the  base  10  are  given  in  parenthesis;  i.e.,  1.00(- 10)  “  1  00  *  10~10 
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Table  2.  H;— 0.  Elementary  Reaction  Equilibrium 
Constants  Jl 


i'b» 

k  -  k,/kr 

- - 

-  A..„T  ;J  exp 

iCq/T)'-'  j 

;  Acq 

| 

I 

K 

- 1 

1  C 

i 

6.131-01) 

i 

-3.991-02) 

i 

9  671+02) 

2 

5  521  +00) 

1 

-3  301-01) 

!  2.861+04)  ' 

3 

3  321+03) 

-6.891-01 1 

1.891  +04)  , 

4 

1  131+02) 

—  5  551—01 ) 

2  '6(  *-04) 

5 

3  461  +02) 

-6  571-01) 

’  091  *-03 1 

6 

1.681  +05) 

-1  30(  +00) 

3  371  +04) 

-T 

7  881+04) 

! 

-1  02  <  -i-  00 ) 

3  831  +04) 

8 

5 091—02 ) 

1.841-01) 

7  801  +03)  j 

4 

1  1 . 671—01 1 

3.601-01) 

9  661  +03)  ! 

10 

3.121-02) 

i 

1.441-01) 

8  *-03) 

11 

-  3.531-01) 

-1  731-01) 

3  63.  +04)  ' 

12 

3.341+01) 

—  5.551—01 ) 

1  481  +04) 

1  13 

2.301+02) 

i 

-3.291-01) 

1.941  +04) 

i  14 

1  781+02) 

i 

-5.111-01) 

2.661  +04) 

1  15 

1  081-02) 

l 

i 

3.771-01) 

2.151  +04)  , 

i  16 

|  1711  +  31) 

-1.621-01) 

!  1.461+04) 

1  17  ; 

2.721-0 3) 

3.991-01) 

1  8.691+03)  ; 

1  18 

9.551+00) 

1 

-2.961-01) 

1  2.751  +04)  I 

;  i9 

4.061-02) 

5.201-01) 

!  4.271+04)  1 

20 

6.761+02) 

-6  401—01 ) 

;  6.101  +03) 

;  21  - 

2.631+02) 

-6.481-01) 

1  4.691+04)  | 

22  ! 

1.801-24) 

-1.711-01) 

>  5.221  *-04)  ; 

'  23 

9  151-26) 

1 

1.301-02) 

6.001+04)  i 

!  24  ! 

3.411-25) 

i 

1.531-01) 

1  2.361+04) 

■  25 

1  361-29) 

8.951-01) 

|  2. 60 1  *-04)  ; 

i  26  - 

2.941-24) 

-1.311-01) 

1  5.121  +041 

27  | 

8.921-28) 

1 

l 

5.571-01) 

i  3.231+04)  i 

I  28 

1.031  —  26) 

2.381-01) 

|  5  991+04)  1 

!  29 

2.661-28) 

6.541-01) 

j  1.281+04)  ; 

J'  Equilibrium  sonsums  jrc  given  in  ihe  non- Arrhenius  lor- 
mat  used  lor  ihe  individual  reaction  rate  constants  and  are 
computed  to  cover  the  .'00* .1000  K  range 
'  V  refers  to  the  reaction  index  <ti  Table  I 
K,  is  given  as  the  concentration  based  equilibrium  wonstant 
1  Kv>  and  ts  dimensionless  except  for  equilibria  involving 
association-dissociation  where  it  is  given  the  anus  m 
cm  /molecule 
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Table  3.  Fuel-Lean  Induction  Times  in  the  high-temperature, 
low-pressure  (HTLP)  Region. 


Test 

No. 

Initial  Mixture 
Composition 
(Moles/Liter)'al 

[Hsl 

[Os] 

[Ar] 

1 

2.991-4) 

9.671-3) 

0 

1200 

51 

50 

2 

2.991-4) 

9.671-3) 

0 

1600 

8 

It 

3 

2.991-4) 

9.671-3) 

0 

2000 

2.5 

4 

1 

4 

7.481-5) 

2.991-4) 

9.601-3) 

1200 

455 

5 

7.481-5) 

2.991-4) 

9.601-3) 

1600 

88 

84 

6 

1  501—4) 

OO 

O' 

9.221-3) 

1600 

44 

43 

|  7 

1.501-4) 

5.981-4) 

9.221-3) 

2000 

14 

16 

iji  Exponentials  to  the  base  10  are  gixen  in  parentheses  le  .  I  OOt-lOl  —  I  00  »  10 

'hi  Esperimenial  data  were  obtained  from  ihe  correlation  function  of  the  incident  shock  data  published  by 
D  R  While  [  54) 

•  c >  Calculated  data  presents  r  as  the  time  at  which  ihe  H  atom  concentration  reached  a  maximum 


Table  4.  Induction  Zone  Exponential  Growth 
Parameters  for  [HO]. 


HsOsAr 

Temperature 

(°K) 

Pressure 

(atm) 

7 

* 

Calc'jl 

Expt,b' 

4:1:95 

1140 

0.62 

11.88 

11.31 

1840 

0.57 

5.56 

6.05 

3:3:94 

1124 

0.62 

10.54 

11.72 

1780 

0.27 

6.00 

4.46 

1:4:95 

1165 

0.64 

9.72 

13.96 

1 

1802 

0.28 

5.72 

6.16 

1:8:91 

1230 

0.72 

8.55 

10.53 

1700 

0.27 

5.58 

5.84 

1  Induction  -lines  -  were  calculated  using  Jechimnwski  and  Houghton  s 
wniermn  d  !HOl  •  l  ')  <  id  mole/ liter  so  determine  :he  end  <>i  he 
induction  nertod 

h,Jjcnimovn»iu  md  Houghton  have  esumuicd  :he  experimental  error  m 
-«  n  ne  .ihout  It)" 
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Table  5.  Comparison  of  Chemical  Equilibrium  Compositions 
obtained  by  Calculation. 


Variable 

NASA"*’ 

Mechanism"” 

X  (Mole  Fraction) 

H, 

4.82(— 02)<c’ 

Os 

1.86(— 02) 

H 

2.38( — 02) 

0 

8.98  ( — 03) 

HO 

2.77  ( — 02) 

HO, 

HsO 

HiO-. 

o3 

2.601-10) 

p(g/cm3) 

1 .43  ( — 05) 

1 .45  ( — 05) 

P(atm) 

1.31  (-01) 

1 .33  ( — 01 ) 

T  (°K) 

2607 

_ 

2634 

ij i  Dale  computed  using  Gordon  end  McBride  modification  of 
N  ^S  \  Program 

' b  •  Data  calculated  using  detailed  mechanism  of  Table  t  to 
generate  time-dependent  ordinary  differential  euuations 
'cl  Read  tabulated  data  in  scientific  notation  as-  I  OOt  —1)1)  - 
I  00  «  10'"' 
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Fig.  1  —  Explosion  limits  for  a  HyOy  2:1  gas  mixture  in  a  KCl-coated  static  flow 
reaction  bulb.  < Adapted  from  Reference  II.) 


LOG  (kr) 
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Fig.  4  —  Comparison  of  calculated  and  experimental  reaction  rate  coefficients  for 
0  FUO  —  HO  +  HO 
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Fig  12  —  Computed  isobaric  contours  of  reaction  induction  time.  -  ,  as  a  function  of  temperature. 

T.  for  HyOy  2:1  gas  mixtures 


(09S rl)  'x 
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d 


.1 


P TOTAL  (atm) 

h*  15  —  Computed  isothermal  contours  of  reaction  induction  time,  r  ,  js  a  function  of  pressure. 
P  ,,r  u  .  lor  a  stoichiometric  H?  -  0;  gas  mixture 
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dlH.O] 

Fia.  18  —  Computed  isothermal  plot  ol  reaction  rate.  - f — •  as  a  lunction  ol  pressure.  P.  above 

dt 

the  2nd  explosion  limit  lor  a  H;:Oy  2:1  eas  mixture 
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V 

Fig.  20  —  Comparison  of  computed  and  experimental  temporal  behavior  of  reaction 
species  concentrations  as  a  function  of  the  reaction  progress  variable,  v 


53 


NUMBER  DENSITY 
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0.158  ATM  H2:02:Ar  /1:3:96 
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0.158  ATM  H2:02:Ar  / 1:3:0 
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